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Abstract 

The vacuum energy density or free energy of a free charged Bose gas at non- 
zero densities is studied in the context of the debate about Multiphcative Anoma- 
lies. Some ("-function regularised calculations of the free energy in the literature 
are reexamined, clarified and extended. A range of apparently distinct answers 
can obtained. Equivalent dimensional regularisation results are also presented for 
comparison. I conclude that operator ordering and normal ordering are not respon- 
sible for these differences. Rather it is an undesirable but unavoidable property 
of C-function regularisation which leads to these different results, making it a bad 
scheme in general. By comparison I show how dimensional regularisation calcula- 
tions give a consistent result without any complications, making this a good scheme 
in this context. 



1 Introduction 



The vacuum energy density in QFT (quantum field theory) is a non-trivial object which 
leads to interesting physical phenomena such as inflation or the Casimir effect. When 
using path integrals to calculate it, one often encounters terms of the form lndet(A^^) 
where A is the propagator for some field. However, in QFT this is the determinant 
of an infinite dimensional matrix, and this is usually infinite. It is therefore extremely 
interesting to find that the "Multiplicative Anomaly" a{A, B) 

a{A,B) := Indet(AS) -lndet(A) -lndet(5), (1.1) 

need not be zero for two pseudo-differential operators A and [1, 2, 3]. The term 
anomaly in this context is used to indicate the failure of a familiar algebraic expression 
rather than a breakdown of a classical symmetry in the quantum theory. Since the 
individual terms in (1.1) are usually naively infinite in QFT, one must regulate before any 
serious discussion. A subscript on a quantity will be used to indicate the regularisation 
scheme used e.g. qr, [lndet(A)]R for some scheme R. The first examples of a non-zero 
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Multiplicative Anomaly in QFT was given by Elizalde, Vanzo and Zerbini [2] using (- 
function regularisation of Dowker and Critchley [4] and two free scalar fields in four or 
more even dimensional Euclidean space time. The Multiplicative Anomaly of (1.1) is 
related to the Wodzicki residue [5, 1] which is important for non-commutative geometry 

[6]} 

There have been two types of extension to the original work of [2]. Firstly, one 
can investigate different renormalisation schemes. In an earlier paper I showed that 
loop momentum cutoffs, including space-time lattices and dimensional regularisation — 
schemes common in particle physics — are all free of this Multiplicative Anomaly [7]. 
Rather, non-zero Multiplicative Anomalies are a feature of the whole family of Schwinger 
proper time regularisation schemes of which (^-function methods are just one example 
[8]. The second direction has been to show that the original Multiplicative Anomaly of 
(1.1) is just one of many algebraic identities which ^-f unction regularisation expressions 
fail to obey [7, 9, 12, 13, 14, 15, 11, 10]. If one is using C-function regularisation, the 
importance of these Multiplicative Anomalies is not in doubt from the point of view 
of mathematical consistency. Furthermore there are several types of problem where ^- 
function regularisation methods are probably the best method to use, for example Casimir 
effect problems or problems of QFT in curved space time. The important debate is 
whether or not the extra terms due to Multiplicative Anomalies lead to new physics. A 
more extensive summary of Multiplicative Anomalies is given in appendix B. 

One of the best examples to study in the discussion of physical meaning of Multiplica- 
tive Anomalies is that of the free Bose gas. This is because even in a free theory with one 
complex relativistic field, the U(l) symmetry can be broken and Bose-Einstein conden- 
sation then takes place [16]. In this model, the calculations of Ehzalde et al [9] (see also 
[10]) of the free energy, using ^-function regularisation methods, gave a result which has 
extra terms when compared to standard results found in texts such as Kapusta [16]. The 
extra terms of [9] were chemical potential but not temperature dependent, and they led to 
alterations of the critical temperature. Furthermore in [9] these extra terms, which come 
from the Multiplicative Anomaly, were shown to be essential for mathematical consis- 
tency. Elizalde ct al. suggested that normal ordering might be relevant in understanding 
why these extra terms appear. McKenzie-Smith and Toms [12, 13] then looked at the free 
Bose gas in (^-function regularisation with particular reference to canonical calculations 
and normal ordering rather than path integral calculations alone. They concluded that 
there is no problem if one follows canonical definitions. 

Unfortunately there are mathematical inconsistencies between some of these (^-function 
regularisation results for the free Bose gas. I will reexamine the claims and calculations 
of Elizalde et al. [9], and McKenzie-Smith and Toms [12, 13] within ^-function regu- 
larisation, and will show how to fit them into a consistent mathematical framework. 
However, I confirm the mathematical results of [9] which show that one can reasonably 
obtain a number of different results for the free energy of the free Bose gas when using 
(^-function regularisation. However I differ from previous authors in my interpretation of 
the physics contained in these results. I conclude that canonical methods, and normal 
ordering in particular, offer no explanation for the genuinely different answers obtained. 
All of the many formal starting points are equally good. Rather my explanation is that 
in certain regularisation schemes, such as ^-function regularisation, the UV divergences 
are controlled by functions whose form depends on physical parameters such as chemical 
potentials. Equivalently the renormalisation scale in ^-function regularisation depends 
on some of the physical parameters of the problem. Hence the comparison of results 

^See Connes [6], page 307. 
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calculated in (^-function regularisation at different values of physical parameters is then 
extremely difficult. I conclude that there is no new physics of the free Bose gas in the 
non-standard results provided by some C-function regularisation calculations. 

However, in order to discuss non-standard results, I must also construct a standard. 
Thus to compare the (^-function regularisation results for the free Bose gas of [9, 12, 13] 
with the free energy results found in the literature such as [16], it is convenient to use a 
standard particle physics regularisation scheme. Dimensional regularisation, as defined 
in the particle physics literature^ [23, 17], is used here for several reasons. Firstly it is 
widely used in practical calculations [17]. Secondly it has no Multiplicative Anomalies 
when used in a sensible and correct manner [7], and so is in the same position as sev- 
eral other standard schemes such as lattice regularisation and simple momentum cutoffs. 
Lastly, dimensional regularisation and ^-function regularisation regulate by inserting a 
non-integer power of a polynomial of loop momentum, so one often finds that the same 
standard integral is required for both schemes. Therefore I will discuss the relation of 
(^-function regularisation calculations of free Bose gas to the standard results obtained 
with dimensional regularisation and in doing so make contact with results found in the 
literature such as are given in [16]. 

In the next section I examine how various formal, and therefore strictly meaningless, 
expressions for the free Bose gas are derived. I give precise well defined mathematical 
expressions only in the following section by carefully defining the regularisations. In 
section 4 I study the relationships between the precise and well defined forms of the free 
energy. Section 5 is devoted to looking at the physical content of the expressions and my 
conclusions are given in the final section. I have tried to put as many technical details 
in the first appendix, but because the confusion in the literature I feel it is important to 
specify the precise mathematical approach used here. Appendix B puts the multiplicative 
anomaly found in the free Bose gas in a wider context. 

2 Formal expressions for the free energy 

The free Bose gas is a heat bath of charged non-interacting scalar particles. The dynamics 
are described by the usual Klein-Gordon Lagrangian for a free relativistic complex scalar 
field=^ 

|a^$|2 _^2|^|2. (2.1) 

For simplicity I work in four-dimensional Euclidean space, of which the spatial dimen- 
sions have volume V, as results such as [2] show that Multiplicative Anomalies are often 
non-trivial in such space-times. The Lagrangian (2.1) is invariant under a global phase 
transformation which is associated with the conservation of particle minus anti-particle 
number^. To describe the statistical average I use the density matrix of an equilibrium 
grand canonical ensemble, for which the temperature and the chemical potential of the 
heat bath (intensive variables) are specified rather than the energy and charge (exten- 
sive variables). The Euclidean approach to thermal field theory encodes this very simply 

^Confusingly, some of the C- function regularisation literature, such as [11], refer to some Schwinger 
proper-time regularisation schemes as "dimensional regularisation" . This is not the same as the scheme 
widely used in particle physics, as noted by Ball [8] and discussed in appendix A. 3 around equation 
(A.20). 

^The condensed matter limit was considered separately in [12]. It is contained in this analysis as the 
TO~/Lt3>'r^TO — /i>0 limit [19] though regularisation prescriptions may differ from that used in [12]. 

"'The number of particles and the number of anti-particles in each mode is also separately conserved 
but these symmetries are lost in interacting theories so I do not consider them here. 



and this is sufficient for our purposes. I will use tfie approacli in wliicli tlie temporal 
direction is made periodic with length f3 while the effects of non-zero chemical potential 
are included by working with an effective Hamiltonian, corresponding to the Lagrangian 
density^ [16, 18, 19] 

= (2.2) 

where 

K± := ((94 ± /x)^ + a;^ (2.3) 
uj^ ■= -v^ + ml (2.4) 

— * 

The V acts only on the three spatial directions. The Euclidean temporal derivatives, ^4, 
are just shifted by the chemical potential fi (a real parameter) corresponding to a real 
shift in the origin of Minkowski energies [19]^. Consider the partition function Z and the 
associated free energy density F 

Z := Tr{e-^(^-'^^)}, F:=--^ln(Z). (2.5) 

I will keep \^\ < m which, in a free theory, means working in the symmetric phase [16]. 
Likewise, I keep T > as at zero temperature all particles will be in the ground state so 
1/^1 = or m. Since the theory is free, one may quickly obtain an expression for F but 
it is crucial to examine the familiar steps, given that the failure of supposedly familiar 
algebraic identities such as (1.1) is at the centre of the Multiplicative Anomaly debate. 

Before I look at these calculations, note immediately that in this section I am writing 
formal expressions, by which I mean they have not yet been regularised, so they are 
naively infinite and therefore strictly meaningless. In doing this I am merely reflecting 
standard QFT procedures and will correct this later sections. This distinction, between 
formal infinite expressions and their meaningful regularised counterparts, is however not 
always clearly maintained in the literature. One of the aims of this paper is to indicate 
clearly when one is performing formal manipulations, and when rigorous manipulations 
of finite objects is being performed. 



2.1 Path integral approach 

When working with the path integral in this model, one would normally exploit the fact 
that the field is complex and work in the one- dimensional complex irreducible represen- 
tation of U(l). The path integral for Z is 

Z ^ J D^D^''exp{- J d^x jC^[^,^^} (2.6) 

'"'The alternative Euclidean method puts both temperature and chemical potential in the boundary 
conditions [18, 19]. 

^The meaning of Hermitian conjugation has been modified in the appropriate manner for Euclidean 
theories. 
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so the usual results for complex Gaussian integration give^ 



F,K+ := -^lndet(ir+) = -^Trln(i^+). (2.7) 

Here the trace Tr with a capital T and lower case r indicates that the trace is over the 
four dimensional Euclidean space so that in coordinate representation 

Tr= r dU f d^x (2.8) 
Jo Jv 

If 1 work with the reducible complex two dimensional representation of U(l), i.e. use 
a vector ^' = ($, $^), then the Lagrangian is given by 

= 1#D-^^ (2.9) 



where 



The path integral for Z is 





K_ 



(2.10) 



j exp{- J (2.11) 

and Gaussian integration gives 

:= -^lnDET(D-^) = -^TRln(D-^) (2.12) 

where TR and DET are now taken over both the two-dimensional field space and the 
infinite dimensions coming from the four-dimensional space-time. One can take the de- 
terminant over the two-dimensional field space [21] to give 

F^A := -^Trln(^) (2.13) 

where 

A^K+.K_. (2.14) 

One often sees this calculation expressed in terms of two real fields, the real and 
imaginary parts of the complex field ^ — {(f)i + i02)/-\/2, e.g. [16]. Writing in a real 

— * 

vector notation = (0i, ^2) gives 

= ^^G-V (2.15) 

^One often sees the Indet notation in C-fmiction regularisation work rather than Trln. However, 
neither has any proper definition unless regularised, and I am free to choose an appropriate definition, 
pace [20]. I will define them to be equal to the same regularised expression. In any case, all practical 
regularised expressions I know involve a sum or integral, i.e. they look like a trace and do not contain 
products which might remind one of a finite determinant. After all, the most familiar definition of 
a Riemann C- function involves a sum over positive integers, rather than the alternative formula of a 
product over prime numbers. I will therefore tend to use the Tr In notation, though I treat the two as 
equivalent. 



(2.16) 



where 

-dl - + - /x^ iidi 

—/j^d^ —84 — + — /x^ 

The path integral for Z is 

Z = j D$exp{- j d^xC^i^l). (2.17) 

In this case the propagator has off diagonal elements proportional to but the eigenvalues 
are the diagonal entries of D in (2.10) and the eigenvectors are $ and These complex 
fields are of course the U(l) charge eigenstates and are therefore the most appropriate 
basis for unbroken symmetry problems. In any case one obtains the same result as before 
(2.13). 

One can now look at a strange variation of the calculation of F^a- One can take (2.17) 
and formally one can factorize it into two integrations 

Z^Z+.Z^, Z± = J D^D^'' exp{- J d^x^^^'K^^} (2.18) 

to give the formula 

Now I see that this seemingly innocent variation with an initial factorization in the 
path integral, is producing a result which is nothing more than a factorization of the 
quartic operator A of the expression F4A (2-13) into two quadratic operators. However, 
it is exactly the success or failure of this sort of factorization which is measured by the 
Multiplicative Anomaly (1.1), i.e. 

a{K+,K_) = 213V {F^K - F^a) . (2.20) 

Thus a non-zero Multiplicative Anomaly says that I can not do these trivial maneuvers 
in the path integral such as described here. Since so much work is based on such formal 
manipulations of the path integral. Multiplicative Anomalies have potentially very grave 
implications. 

Returning to the formal expressions for the free Bose gas free energy, note that there 
are many other factorizations of the operator A. I will follow [9] and consider just one 
other, namely 

L± := -^l + {uJ±^if (2.21) 

where uj was given in (2.4). I could therefore equally well define the free energy density 
as 

Again, the Multiplicative Anomaly expresses the failure of algebraic identities needed to 
relate F^a to F^i and we have from (1.1) 

a(L+,L_) = 2f3V{F4L-F4A) (2.23) 

This shows again that the Multiplicative Anomalies encode differences between formally 
equivalent definitions of the free energy of a free Bose gas. 

Several other combinations also exist but I will focus on just the L and K factoriza- 
tions. One common feature of all these expressions is that they involve four dimensional 
traces, equivalently are evaluated using integrals over energy and momentum, and that 
they are most naturally obtained within a path integral approach. 



2.2 Canonical approach 



The free Bose gas gas is also easily obtained using canonical methods and, as emphasized 
by McKenzie-Smith and Toms [12, 13], this leads to additional familiar forms. 1 will 
start with the Hamiltonian and charge operators for a relativistic free Bose gas. Unlike 
the path integral with its c-numbered fields, I must specify an operator ordering. My 
first definition of a Hamiltonian and charge operator are ones which have been normal 
ordered in the conventional zero temperature sense^, i.e. annihilation operators to the 
right, which I will denote with a subscript A^, namely 



Hn y (2^(^.4% + ^^4^^)' (2-24) 



with dispersion relation uj^. Right at the start note that 1 am working on mass shell 
with integrals over three momenta, not traces over four space-time coordinates. Thus 
the expressions obtained will have a trace over the three dimensional spatial coordinates, 
denoted with lower case letters, as 

tr= / d?x, (2.26) 



V 



to distinguish it from the full four-dimensional trace Tr of (2.8). 

The vacuum state of the Fock space associated with the and 6g operators is 

%|0)='^dO) = (2-27) 

This vacuum state will be the true physical vacuum provided we avoid symmetry breaking, 
i.e. if I choose > and if I have low charge densities so that |//| < m so there is no 

Bose-Einstein condensation. 

With these definitions, I see that the expectation values of both and Qat in the 
vacuum state are zero, 

(0|#jv|0) = (OlQjviO) = 0. (2.28) 
The density matrix in terms of these normal ordered operators is just 



:= exp{-/3(i/w-/igiv)} = exp{-/3 / y ((o;^ - /i)aL% + {uj^ + /x)6t6g 



(2.29) 

Another ordering is the symmetric one, which I will denote with a subscript 5", namely 

^5 / ^ (^^^(4% + %4) + t'S'h + ' (2.30) 

- \^^k^WS\ ■ (2-31) 



Such a symmetric operator ordering is often the ordering found necessary to match path 
integral results, though a path integral is written in terms of commuting objects. Note 



^Interestingly, at T > or /i 7^ the most appropriate normal ordering is different and depends on T 
and /z explicitly [22]. 



that, at least formally, the expectation value of this Hamiltonian in the vacuum state is 
now not zero, but is equal to a half a;^ per mode in the system. On the other hand the 
expectation value of this charge operator is still formally zero, 



N 



|0>=2(/ 



(2^3 2 



1 



(OIQatIO) -0. 



(2.32) 



The symmetric ordering density matrix is then just 



— exp{—(3 



{2ny 



1. 

2K 



/^)(4«fe + «fe4) 



(2.33) 

^K + /^)(^P^+¥p)(2-34) 



One can imagine other operator orderings. In particular, the only way to get charge 
operators with a non-zero vacuum expectation value is to change the ordering from that 
of Qn by different amounts for the and 6^ operators to give a charge operator equal to 
Qas = Qn + o: where a is a c-number. This is equivalent to changing the definition of zero 
of charge. It will not effect the physics as all such Qas operators still commute with the 
Hamiltonian, but the labels I give to a given physical state will depend on such decisions. 
Each (anti-)particle still carries charge +1 (—1) but now the state with equal numbers 
of particles and anti-particles has charge a. This does not seem useful physically and 
further I shall show that normal ordering has nothing to do with the anomalous terms 
of interest here. It is therefore sufficient to focus on just the normal ordered (2.24) and 
symmetric (2.30) versions. 

The canonical calculation runs as follows. Since it is a free theory, the Hilbert space 
is a direct product of the single mode states, that is a direct product of the Fock spaces 
associated with the single oscillators %(A;) or bf:{k) for all possible k. In practice this 
means that for the symmetric density matrix (2.33) 



= Tr{e 



-PiHs-lxQs) 



} 



(2.35) 



n t^a,fe{exp{-/3^(^fe - Ai)4%}) 



X 



ntv{cxp{-/3^(a;,^ + /x)&t6,^}} 



(2.36) 



where here the traces, tr, are over the appropriate single oscillator Fock spaces. The 
direct product of the Hilbert space has resulted in a factorisation of terms in this formal 
expression. Standard calculations then give 



3/i 



E 

± 



\H^k ± + ^tr{ln[l - exp(-/3(a;g ± 



(2.37) 



In exactly the same way, the normal ordered form for the density matrix (2.29) gives a 
second formal expression 



SAT 



E 

± 



-tr{ln[l-exp(-/3(^,^± /.))]} 



(2.38) 



o 



with the (3 independent zero-point energy term missing as one would expect. 

However, there is yet another formal expression discussed in the literature. This I will 
call F3 where 

F3 := tr{a;J + ^trln[l - exp(-/3(a;^ - //))] + ^trln[l - exp(-/3(a;g + /i))](2.39) 

This has ± A* factors in temperature dependent terms but not in the temperature inde- 
pendent terms, so is perhaps not the most natural form to arrive at by direct calculation 
though physically it is quite intuitive. The obvious way to obtain this third form F3 (2.39) 
is by making algebraic manipulations of the formal form (2.37). Later, inspired by the 
work of [13], I will find different answers in ^-function regularisation for these three forms, 
(2.37), (2.38) and (2.39). Normal ordering gives a genuine formal difference between the 
formal expressions (2.37) and F^n (2.38) but not between (2.37) and F3 (2.39). 
The difference between ^3^ (2.37) and F3 (2.39) is a matter of the Multiplicative Anomaly 

ashift(^, a) ■= 2tr{A} - tr{A + a} - tr{A - a} (2.40) 

This expresses the failure of two basic algebraic identities (B.5) and (B.6) discussed in 
appendix B. In the case of the three dimensional expressions F3 and ^3^, the difference 
is related to the Multiplicative Anomaly (2.40) through 

F3 = + ^ashift(a;g, 1^) (2.41) 

Thus the unexpected differences between various C-function regularisation versions of 
the formally UV divergent canonical calculations F3 and F^^, as found in the literature and 
to be discussed below, are not caused by normal ordering but failures of basic algebraic 
identities. It is the difference between F^n normal ordered UV finite expression and F^^ 
UV infinite expressions which is a matter of normal ordering and the zero point energy 
in this model. In this I differ from both [9] and [12, 13] who emphasise the role of normal 
ordering and canonical methods in their resolutions for the problems. 

3 Free Bose gas results 

I will now turn from infinite formal expressions to finite regularised ones. Only with 
the latter can one compare the results for physical quantities obtained using different 
calculational schemes. I will denote regularised quantities by adding a further subscript: a 
C to indicated a (^-regularised expression, or a e to indicate that dimensional regularisation 
was used, e.g. F4K( is the ^-function regularised form of F^k- 

Before looking at the results in detail, I will make a few comments about the way 
I approached the four-dimensional finite temperature calculations of (2.13), (2.19) and 
(2.22) where there are energy variables to be summed over. As noted above, I use a 
Euclidean approach to thermal field theory where the energies are the discrete Matsubara 
energies but the chemical potential is encoded directly in the Lagrangian, propagators 
etc. rather than in the boundary conditions [18, 19]. There are standard methods for 
performing the Euclidean energy sums, for instance using contour integration methods. 
In this case the integrands are logarithms or non-integer powers of polynomials of energy 
so more care is needed than with simple Green functions with their integrals of rationals 
of polynomials. Nonetheless the discussion in the appendix A.l shows that, as with Green 



n 



functions, one can separate any calculation into two pieces 



/ + 00 
dhg{h). 
-oo 



(3.1) 



(3.2) 



where in (3.1) g contains all the /i dependence. In (3.2) Fp contains all the explicit 
temperature dependence in factors of exp{/3(a; ± /j)}^. The Fp are also guaranteed to 
be UV finite. The first piece, Fq, has no explicit temperature dependence and is just 
the original expression with the energy sum replaced by a Euclidean energy integral 
from minus infinity to plus infinity. However, the zero temperature piece Fq will be 
chemical potential dependent. As the problem of Multiplicative Anomahes is all about 
the behaviour of UV divergences, the focus will be on the temperature independent but 
possibly n dependent UV divergent parts Fo(a*)- The major question is how to implement 
the regular isat ion of the UV divergences in these temperature independent, chemical 
potential dependent terms. 

Various terms crop up again and again so I will define some useful functions. 



/_i(m) /_! 
/o(m,M) :=/o 

Y{m,i^) := Y 

X{m,n) := X 



m 



327r2' 





In 
















2/-i(m) 







m 

32^ 



In — 



m 
M2 



1^ 
3 



(3.3) 
(5-4) 
(3.5) 
(3.6) 
(3.7) 



3.1 4D (^-function regularisation results. 

The details about how to implement (^-function regularisation in this simple case are 
given in appendix A. 3. 1 have chosen to write these (^-function regularised expressions so 
as to make clear the close relationship to the dimensional regularisation results of later 
subsections. This will simplify the comparisons 1 make in section 4 but it means that 
my ^-function regularisation expressions are not always exactly the same as those in the 
literature, but the relationship is trivial (sec appendix A. 3). My ^-function regularisation 
versions of the formal expressions (2.13), (2.19) and (2.22) are 



MAC 



{kj+u^l + l^r-Au^y 



± 



d^k 



M4 

{k^ ± ifi f + Lol 
M2 



-s/2 



(3.8) 



(3.9) 



^For Bose or Fermi statistics and m, > such arc always zero at zero temperature. However the 
charge density is the physical parameter, not /U, and this makes an implicit function of temperature. 
See appendix A.l for further comments. 



where I have chosen to ensure that in all expressions the regularisation scale always 
appears as M^* and the first term in a small s expansion is always 0{l/s) to maintain 
a close analogy with dimensional regularisation's M^^ and 0(l/e). Using standard tricks 
it is straightforward to do these integrals and I find for (3.9) and (3.10) 

F4i^C - /-i(m)i + /o(m,M) + /^(m,/x,/5) (3.11) 
s 

F4LC = U(m)- + fo(m,M) + 2Y(m,ii) + fp(m,ii,(3) (3.12) 
s 

Note that both of these are calculated directly from the expressions quadratic in energy- 
momentum used to define them in (2.19) and (2.22), that is no use of Multiphcative 
Anomalies was made in their calculation. The surprise is that ^4^,^ is different from F^kc 
by just the sort of (3 independent, /j, dependent term I will be talking about in the context 
of Multiplicative Anomalies. The expressions are essentially the same^° as given in [9]. 

Now one can consider the expression ^4/1^ based on quartic energy-momentum terms. 
The calculation has been done in three ways. First, one can work directly with the quartic 
terms and [9] obtain 

F4A( = /_i(m)- + /o(m,M) + r(m,//) + ^(m, //,/?) (3.13) 
s 

In the second and third approaches to F^ac, I start from the definition of F^ac in terms 
of Trln over quadratic operators and an additional anomaly term, as given in (2.20) and 
(2.23). Results for the Trln over quadratic operators were given above in (3.11) and 
(3.12). However, one sees from the definition of the Multiplicative Anomaly (1.1) that 
to get the result for F4^ac these cases 1 must calculate the relevant Multiplicative 
Anomalies, as indicated in (2.20) and (2.23). These types of Multiplicative Anomaly can 
be calculated directly from a formula based on the Wodzicki residue [5] of the theory of 
elliptic pseudo-differential operators. Details are given in appendix A.4. In any case, I 
find that 

a{K+, K_) = -a{L+, L_) = -2pVY{m, /x) (3.14) 

which agrees with earlier results of Elizalde et al. [9]. It is then clear that all three methods 
of calculating F4AC give the same answer and therefore I disagree with the suggestion of 
an inconsistency made in [13]. 

The fact that all three approaches give the same answer is fundamental. It is a matter 
of mathematical definition that all three methods — direct, via K± and its Multiplicative 
Anomaly, and finally via L± and its Multiplicative Anomaly — are calculating the same 
object. However, if we had forgotten the Multiplicative Anomaly then I would have an 
inconsistency in our results for ^4^^. In this sense the existence of the Multiplicative 
Anomaly is not a problem but is absolutely essential for mathematical consistency.^^ 

Note that while there is complete consistency in the F^ac results, the Y factor comes 
up in different ways in different approaches. When using F^kc (2.20) approach the Y 

^'^The only differences arc the sign in front of the /o term in (3.11) (equation (44) of [9]) and a factor 
of two in (3.12) (equation (38) of [9]). 

^^Of course, if you accept that Anomahes exist then this shows how to calculate them without knowing 
the Wodzicki formula. Calculate -F4AC, and F4KC directly as done in [9], and then one can deduce the 
result for a{K+, K-). 
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comes solely from the Multiplicative Anomaly a{K^, K^). When starting from the L 
based form ^4^,^ of (2.23), there is a F term intrinsic to F^li^ and another one from the 
Multiplicative Anomaly a(L+, L_) which contribute to the final Y factor in F4a<: (3.13). 

Of course, while mathematical consistency of the ^-function regularisation may require 
non-zero Multiplicative Anomalies, this appears to cause major problems to the physical 
interpretation. I will return to this later when all relevant results have been acquired. 



3.2 3d (^-function regularisation results. 

The details of the finite temperature aspects and (^-function regularisation are as for the 
four- dimensional calculations above (also see appendices A.l and A. 3). Thus I define the 
^-function regularisation forms as follows 

:= ^E/(05(-±'')""'>^" (3-16) 

and Fsiv is UV finite so needs no regularisation. I choose s so that the renormalisation 
scale appears to the same power as in the four-dimensional calculations. Direct calculation 
then gives 

FsN = U (3-17) 
Fsc = /-i^ + /_i(ln(^^+V'(3)-V'(-l/2))+/^ (3.18) 

= /-i^ + /o + /-i(V^(l)-^(-l/2)) + /^ (3.19) 
FsK = F3c + 2r(m,/x) (3.20) 

Thus I see that in ^-function regularisation 
ashift,c(^, 1^) = 2tr(c^) - tr(cu - /x) - tr{u + ji) = 2 (Fg^ - Fg^c) = -4:pVY{m, /i). (3.21) 

The Wodzicki residue formula (A. 24) is only useful for TV In expressions. Attempts to use 

it here forces one to try exponentials of well behaved operators in the formula (A. 24) but 
these are not suitable. However, I am able to get the form of this Multiplicative Anomaly 
directly, if not the overall factor, by using a conjectured generalisation of the Wodzicki 
residue formula (A. 30). 



3.3 4D dimensional regularisation results. 

It is useful to compare the ^-function regularisation results against results of a method 
more often used in particle physics, so I will use dimensional reduction. When using 
dimensional regularisation at non-zero temperature it is important to note that the reg- 
ularisation takes place solely in the spatial integration, as discussed in the appendix A. 2, 
i.e. 

4'-'''k = (27r)2^-^-^ Yl j ^'"''^ (3-22) 
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Let us start with the quartic expression (2.13). The appropriate dimensionally regularised 
form at T > is 



/^^--^M-ln(^) (3.23) 



The key point here is that this expression is finite while e 7^ so I can manipulate the 
integrand and swap the order of summing integrands with integration, for instance 



4Ke 



Thus it is immediately obvious that there are no Multiplicative Anomalies in dimensional 
regularisation as noted in [7]. Denoting all calculations in dimensional regularisation 
based on any of these four-dimensional forms as F^^ I find 

F4e = /-i(m)^ + /_i(m)(^ln(^^^+ln(47r)+V(3)^+/;3(m,/x,/3) (3.25) 

= /_i(m)- + /_i(m)(^ln(^^^j+V^(3)-V^(l)j+/^(m,/x,/3) (3.26) 

- /_i(m)i + /o(m,(47re^)^/2^) + /^(m,/x,/5) (3.27) 

Before moving on, note that is possible to produce failures of algebraic identities in 
dimensional regularisation if one fails to implement dimensional reduction correctly. For 
instance this can be done by regulating the energy sums rather than the three-momentum 
integrals, see appendix A. 2 for further comments. 

3.4 3d dimensional regularisation results. 

Just as I noted for the four-dimensional calculations using dimensional regularisation, 
if one follows the standard application of the dimensional reduction scheme to one of 
the three dimensional forms (A. 7) one quickly sees that there is no problem with simple 
algebraic identities such as (B.5) and hence the Multiphcative Anomaly Oghift of (2.40) is 
zero in dimensional reduction. Thus there is no difference between the two dimensionally 
regularised expressions based on the symmetric ordering, of (2.37) and F3 and (2.39). 
Likewise, the third expression, F^^e is just the UV finite T > term present in the all 
the other expressions encountered in dimensional reduction or ^-function regularisation 
calculations. The dimensionally regulated form of both (2.39) and (2.37) is then 

Fs^e^Fs, := J f-^'k M^' cuik), F^Ne ^ (3.28) 

Prom this I find 

F,,{m,M,,l^,T,e) = /_i(m)^ + /_i(m) (^In (^^^ + ln(47r2) + ^(3)^ + /;{3.29) 

We see that the three- and four-dimensional results (except for the UV finite F^^^) are 
identical with the same choice of e and renormalisation scale Mg. Prom this I find 

Fse{m,M„ii,T,e) = Fi,{m, M„ ii,T,e) (3.30) 



This is not entirely trivial but comes as a result of the careful implementation of the 
dimensionally regularised scheme, e.g. before integration free energies are always propor- 
tional to M^^. This ensures identities such as (A. 8) are always satisfied which helps hnk 
three- and four-dimensional results. One can also make this connection directly. First, 
remember that F^Ae is regulated only in the spatial integration, as shown in (3.23). Do- 
ing the Euclidean energy sum is thus straightforward using the usual contour integration 
methods as all the regularisation is in the spatial momentum. This shows that there are 
no T — iJ, dependent terms, just a T = UV divergent piece and the UV finite T > 
contribution fp}"^ More relevant is that the T = piece of is then a three-dimensional 
integral and is then clearly of the form (3.28), the three-dimensional expression ^3^. 

4 Comparison of results 

I have focused on six different expressions using two different regularisation schemes all 
purportedly for the same quantity - the free energy of a free charged Bose gas. As always 
in QFT, comparison of results is complicated by the need to make sure that in two 
expressions for the same quantity that the physical parameters, upon which the answer 
depends, are defined in the same manner in all cases. Since it is a free field theory, 
there are no obvious corrections to the bare mass which is the physical mass, despite the 
quantum and statistical fluctuations encoded in the calculations. However one also has 
the quartic UV divergences in the free energy density or equivalently in the zero point 
energy density. Thus to discuss the physics in these results one must not only discuss the 
regularisation of the UV divergences, but also the removal of these divergences as part 
of the renormalisation procedure^^. I will postpone the discussion of the renormalised 
results and the physics to the following section and in this section I will compare the full 
regulated results. 

The regulated but unrenormalised (^-function regularisation and dimensional reduction 
results can be compared directly by setting e = s as by inspection one see that all 
the singular terms are then identical. This is not a miracle but comes from the close 
relationship between dimensional reduction and ^-function regularisation and from the 
way I have chosen various optional factors in such a way as to ensure this behaviour, see 
the appendices for further details^^. Also note that the expressions for the Multiplicative 
Anomalies are themselves naively finite and in principle require no renormalisation to 
obtain finite answers, though one needs to implement regularisation carefully to ensure 
this is seen in the final answers. 

However, regularisation alone introduces a new scale, M the renormalisation scale. 
This is not a trivial object. A good practical example of the importance of M is the way 
that in particle physics, a great deal of effort goes into setting the scale for lattice and 
MS dimensional regularisation results for the same quantities [17]. 

The conclusion is that in order to be able to compare the different results for the free 
energy obtained above I must not assume that the renormalisation scales M appearing 
in different calculations are the same. 

^^The T = part of F4K( can be rewritten using the identity (A. 8) as the usual four- dimensional zero 
point energy result for a complex scalar field in dimensional regularisation. 

^^The latter process is needed when extracting the physics in QFT whether or not there are UV 
divergences. 

^^These extra factors, functions g in (A. 7) or (A. 12) in the appendices, can easily be included provided 
one modifies the relationship to s = e/g{e) or similar. Fixing g = 1 still leaves the freedom to manipulate 
the renormalisation scale M. Thus I choose to work with g = 1 merely to remove irrelevant complications. 

1 A 



4.1 Discussion of dimensional regularised results. 

As I will consider these to represent the 'standard' results, I will discuss these first. 
Several aspects are well known but, in view of the complications associated with the ^- 
function regularisation results, I will repeat them and emphasize several aspects for later 
comparison with the ^-f unction regularisation results. 

The canonical expression for the normal ordered expression F3JV leads to the non-zero 
temperature term fp. This is UV finite and hence independent of regularisation. is also 
zero at zero temperature. This accords with the usual idea that normal ordering removes 
the zero point energy term associated with expectation values in the pure vacuum state. 
Thus f/3 appears to be the T > UV finite correction to the zero point energy, and indeed 
we will find this borne out below in other non-normal ordered calculations. The finite 
temperature literature usually focuses on this term alone and leaves the temperature 
independent UV divergent zero point energy terms to one side e.g. [16]. 

The other five expressions for the free energy density, which are all UV divergent, lead 
to identical results in dimensional regularisation. Specifically the four-dimensional results 
-^4Ae(3.23), F4Ke (3.24) and ^4^^ lead to (3.25). The two three-dimensional expressions F^^ 
and Fs^e of (3.28) and (2.37) respectively lead to the single result (3.29), which is identical 
to the four-dimensional dimensional reduction result (3.25) if the same renormalisation 
scale M is used^^. Since equations such as (2.20) show that any differences in results are 
related to non-zero anomalies, my conclusion is that there are there are no Multiplicative 
Anomalies in dimensional regularisation of the free Bose gas model. This agrees with 
the suggestion made in the introduction, based on the criteria set out in [7] , namely that 
conventional dimensional reduction is always Multiplicative Anomaly free. 

The result for these UV divergent expressions comes in two parts. First, the difference 
between T > and T = values (given that < m) is always the same term as in F^e- 
The remaining part, r~^/_i -|-/o (where r is either s or e) is independent of temperature, 
and contains both finite contributions and UV divergences. An important point is that 
in dimensional regularisation this T = part is also independent of chemical potential. 
Since the normal ordered form F^n should have removed the zero point energy, (^3 — ^3^) 
is the zero point energy in dimensional reduction. The zero point energy in dimensional 
reduction is made up of the temperature and chemical independent terms (infinite and 
finite) as one would expect from its definition in terms of a vacuum energy expectation 
value^^. Thus normal ordering merely removes the UV divergent zero point energy, it 
is not the source of any Multiplicative Anomaly. Indeed there are no Multiplicative 
Anomalies in dimensional reduction [7], as discussed earlier, while normal ordering has 
the expected effect in dimensional reduction, clearly showing there is no link between 
Multiplicative Anomalies and normal ordering. If any final proof of this is needed, recall 
that the appropriate normal ordering at T > 0, /x 7^ is not the usual one considered at 
zero temperature [22]. 

One last comment upon the dimensional regularisation results is to note that all of 
the UV divergent forms in dimensional regularisation are equal when the same renormal- 
isation scale. Me, is used. This is a direct consequence of the consistency demanded for 
dimensional regularisation integrals [23] and in particular comes from the result (A. 8) 
which can be used to link the three and four-dimensional results directly with the same 
scale Me. 

^^The fact that the three- and four-dimensional dimensional regularisation results are equal with the 
same renormalisation scale follows from identities in dimensional regularisation such as (A. 8). 
^^One may easily work at T = and 7^ to confirm this directly, as appendix A.l shows. 
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4.2 Discussion of (^-function regularisation results 

The ^-function regularisation produces the same result as dimensional reduction for the 
UV finite normal ordered expression F^n/^. After this, things become more complicated. 

Let us start with the four dimensional form F4K. The C-f unction regularisation result 
for F4K has no terms which depend on chemical potential other than the UV finite 
e.g. no Y factors. This is the same behaviour as all the dimensional regularisation results. 
One can see this result quickly when using the Schwinger trick 

1 1 r°° 

dt f-\- 

a' r(s) io 

on (3.9). There the T = term has an energy integral which becomes a Gaussian with 
a peak at // so a simple shift removes all n dependence from the T — terms. By 
inspecting the results for F4KC it is therefore straight forward to make a link between the 
dimensional regularisation results and this (^-function regularisation result. From (3.25) 
equality is obtained by choosing slightly different renomalisation scales for this ^-function 
regularisation calculation and the dimensional reduction calculations 

F,{m,M„n,T,e^s) = F^Kcim, M^kc, l^,T, s ^ e), (4.2) 
{Ane'^y/^M, = M^kc (4-3) 

The problem is with the other two four-dimensional (^-function regularisation results, 
-^4Ac of (3.13) and F^lc of (3.12). Both differ from F^kc by factors of Y, which is both 
mass and chemical potential dependent and which can not be absorbed as constant shifts 
to the regularisation scale. In fact I see that to obtain equality I must set 



FiAcim, M^Ac fJ^,T,e = s) = F^Kcijn, M^kc fJ.,T,s = e), (4.4) 

e^M4Ac = M4A-C (4.5) 

F^Lcim, M^LC fJ.,T,e = s) = F4Kc{m, M^kc, fJ^,T,s = e), (4.6) 

e2^M4Lc - M4KC (4.7) 



where X{m,n) is defined in (3.6). Thus as < < m I have 1 > e"^ > 0.513 and 
1 > > 0.264, sizable shifts in the renormalisation scales if I want to demand equality 
of free energies in these different C-f unction regularisation schemes. 

Having dealt with the four-dimensional forms, let me now turn to the three dimen- 
sional results using (^-function regularisation (3.18) and (3.20). Comparing with previous 
(■-function regularised results for the four-dimensional calculations, I see that 

F3^{m,M3^,li,T,s) = F4Kdm,M4c,li,T,s), (4.8) 
F3i,cim,M3c,H,T,s) = F^idm, M^c H,T, s) (4.9) 

provided I shift the renormalisation scale for these three-dimensional (-function regular- 
isation results as compared with the four-dimensional scales as 

= M|^exp{V'(-l/2) - V'(l)}, => Mac = ^M^^. (4.10) 

In fact one can derive these relationships directly from the four-dimensional forms, and 
in doing so prove explicitly the relationships (4.8) and (4.9) between the three- and four- 
dimensional ^-function regularisation expressions. This can be done by starting with 
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the four-dimensional form, rewriting the integrand using the Schwinger trick (4.1), then 
doing the energy integration, and finally inverting the Schwinger trick (4.1). 

However one tackles the problem, the results (4.8) and (4.9) show that the three- 
dimensional C-function regularisation forms, as used in [12, 13], hold no new lessons as 
compared to the four-dimcnsional case, which were the focus in [9]. As in dimensional 
reduction, normal ordering gives a difference between normal ordered F^j^^^ and sym- 
metric ordered three-dimensional ^-function regularisation expressions, specifically 
it removes the UV divergences. However the Multiplicative Anomaly problem comes in 
comparing and F^^^, both based on the same symmetric ordering of operators and 
so any difference is not due to operator ordering. Since the three- and four-dimensional 
results can be linked by (4.8) and (4.9), the conclusion is that operator ordering can not 
be responsible for the differences in UV divergent four- dimensional expressions either. 
The slightly more complicated relationship between the renormalisation scales in three- 
and four-dimensional calculations is to be expected in general regularisation schemes, and 
was only avoided in dimensional reduction because of identities such as (A. 8). 

Still, in (^-function regularisation it appears that one can obtain a single result for the 
UV divergent expressions only if we allow a rescaling of the renormalisation scale that 
depends on both mass and chemical potential. By comparison in dimensional reduction 
simple equality of all UV divergent results was achieved. 

Having shown how one can link all these different results, I will now turn to the most 
important question, namely are there any physical differences encoded by these results. 



5 The physics of the free Bose gas 

The results of the previous section are summarised in table 1. 



Quantity 


Coeff. of 

r-7-i + /o 
(T = // ind.) 


Coeff. of 
Y 

(T = /X dep.) 


Coeff. of 

U 
(T > 0) 


Relative 
scale used 


Generic F 


a 


b 


c 


z 


FsNe and F^^c 
All other 

F4KC 

FiLC, 
FiAC, 

F^c 

^3K 



1 
1 
1 
1 
1 
1 





2 
1 

2 


1 
1 
1 
1 

1 
1 
1 


(47re^)-V2 
1 

e/2 
ei-2^/2 



Table 1: Table of the different terms appearing in each expression for the free energy. 
The necessary rescalings of the renormalisation scale relative to the F^kc, calculation scale 
are given. The generic example is therefore F = a(r^^/_i -|- /o) + bY + cfi^ using a scale 
M — zM^KC- r — s or e as appropriate. 



However, before any discussion of physics in these results can take place, one must renor- 
malise the free energy density, F. I will consider two distinct methods of removing its 
quart ic divergences, first a physical subtraction and then an unphysical minimal subtrac- 
tion. 
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5.1 Physical Subtraction 



Consider a free energy difference, AF, defined with respect to a reference temperature 
Tr and cliemical potential /ir, all other physical parameters held fixed 

AF(m, M, T, ij) = F(m, M, T, //) - F(m, M, Tr, hr). (5.1) 

Such a subtraction will not alter the thermodynamics and it is common to work with 
energy differences in many problems. One quickly finds that there are three distinct 
results for AF. For all the dimensional regularisation calculations, and for the ^-function 
regularisation calculations based on F^^ki^i Fsnc -^ac -'- have 

AF, AF^KC = AFsATc = AFg^ = /^(m, T, /x) - /^(m, Tr, iir) (5.2) 

where fjj of (3.7) is the usual T > term encountered in the free Bose gas model e.g. in 
[16]. The (^-function regularisation result based on the quartic operator, ^4^^, picks up 
an extra Y difference term 

AF4AC = AF, + y(m, ^i) - y(m, ^JiR), (5.3) 

while the L factorisation, Fi^^, and one of the three-dimensional results, Fs^^, in Q- 
function regularisation give 

AF4LC = AFa^c = AF, + 2Y{m, //) - 2Y{m, ixr). (5.4) 

These are all UV finite, though in most interacting models a single subtraction would 
not normally leave a finite result in this way. They are also all independent of the 
renormalisation scale M, provided we use the same value for M at T and n as at the 
reference point Tr, fiR. One can check that the physical quantities, various derivatives 
of AF, as a function of the observables, T and the charge density Q/V rather than a 
function of T and /x, do differ depending which result from (5.2), (5.3) and (5.4) we take. 
They are though independent of the reference point Tr, hr. At the same time, all the 
different definitions of F given in section (2) seem to be equally good. Thus there appears 
to be a serious problem in identifying the physics even in this simple free Bose gas model. 
Before suggesting a resolution, let me now consider a different renormalisation scheme. 



5.2 Minimal subtraction 

In minimal subtraction, one just drops the divergent terms. This is a common scheme in 
dimensional reduction and it is also performed implicitly in ^-function regularisation as 
explained in appendix A. 3. Now I find four distinct results, which I denote as F^^. The 
UV finite normal ordered results are unchanged by this renormalisation 

F^l = F^\ = Mrn,T,^i). (5.5) 

Otherwise, the results fall into three types depending on the calculational scheme, just 
as they did for the physical subtraction renormalisation used above. Thus 

i^4x' = /o(m,M) + /^(m,T,/x) (5.6) 
i^4LC = fo{m,M) + f^{m,T,fx) + Yim,fi) (5.7) 
F^C = Mm,M)+Mm,T,i^) + 2Y{m,^) (5.8) 
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All these renormalised free energies based on UV divergent expressions are finite again 
but now all have a factor of /o (3.4) unlike in the case of the physical subtraction results. 
Thus these all depend on the renormalisation scale M. However the /o factors do not 
alter the thermodynamics, assuming the same M is used at all T and /i, so these minimal 
subtraction calculations give the same physical results as when the physical subtraction 
renormalisation is used. The physical results appear to be independent of the renormal- 
isation scheme yet depend on the details of the definition and regularisation of the free 
energy. 

5.3 Resolution of the paradox 

The clue comes from the minimal subtraction results which retain some explicit M de- 
pendence even though this does not effect the physical thermodynamics. It is a reminder 
that when comparing results at different T and /x (or even different m), we are naturally 
inclined to hold this renormalisation scale constant. This is implicit in the physical sub- 
traction calculations (5.2), (5.3) and (5.4) where M is assumed to be the same at T, /x as 
at the reference point Tn, hr. 

But why should the renormalisation scale M be held constant? Well it first appears 
when we cutoff the UV modes at some scale A using a regulating function R{k). In the 
language of this paper, A ~ M/r where r = e in dimensional reduction and r = s in 
^-function regularisation. Even though the interesting physics will be happening at scales 
much less than A and will be encoded in 0(r°) terms of a small r series, there will also 
be a memory of the UV modes in the 0(r°) terms and these will be M dependent. 

Now if we compare a regularised free energy expression for different physical parameter 
values, e.g. calculate AF, it is absolutely essential that we have subtracted these UV 
modes in exactly the same manner in both expressions. If not each expression will have a 
slightly different remnant of the UV modes and the difference in the way the UV sector 
was treated will appear in the physical results and may be finite. It is though a purely 
mathematical artifact, a sign that we did not deal with the UV sector carefully enough 
and that we are not comparing like with like. 

For instance, the simplest way to cut off the UV modes is to use a straight cutoff, 
inserting a regulating function Rx = 9{M/X — \k\) (A — > 0) into all loop integrals^''. This 
will give terms such as 

Fx = I (0I^a(A;) ln(^^ + m') = cM'^ + ... (5.9) 

where c is some number. The UV cutoff scale is at A = M/A^/^ where A is the small 
parameter playing the same role as e in dimensional reduction and s in (^-function regular- 
isation. There is then no doubt that we are removing the UV modes in the same manner 
for all calculations and, provided that the physics was occurring at scales much less than 
A = M/A^/^, no physics would be effected. If we were then to compare expressions for 
the free energy at different temperature or chemical potential in the free Bose gas model, 
e.g. calculate AF, the expected result of [16] would appear, namely it would depend only 
on ffs factors and no cAf^A^^ type term would survive, no extra T oi fi dependent terms. 

Suppose though that one chose the cutoff A(/i) to have some slight dependence on the 
chemical potential, e.g. A*(//) = M^/A^-|-//^ (or indeed any other physical parameter). If 

^''At non-zero temperature, the best approach is to cutoff only the three-momentum, but such detail 
is not relevant to the discussion here. 
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one compared the same regulated expression at two slightly different /i values (e.g. looking 
at a /X derivative of the free energy needed to calculate the physical charge density) , terms 
depending on A would not exactly cancel. For instance making a physical subtraction 
we would find AFa ~ c(//^ — + ... for the simple example mentioned above. A finite 
remnant of the UV mode contribution is left, it is chemical potential dependent and it 
therefore alters the thermodynamics. 

Thus in the case of a simple cutoff, the regularisation procedure is obvious. The cutoff 
A must be taken to be independent of physical parameters. 

In the case of a simple cutoff, making the cutoff A, or equivalently M, depend on 
physical parameters is clearly a bad idea but it is easily noted and avoided. However, 
what if the form of the cutoff function R was /j, dependent rather than the cutoff scale 
itself? For instance 

Rx = 2[exp{(|fc| - M/X)/w} + 1]-^ (5.10) 

is a theta function with its sharp jump at M/A smoothed over a region of size w. If 
the width w varies slightly with fj, then this will generate spurious /i dependent terms 
from high energy k ~ M/A modes when comparing free energy calculations at different 
chemical potentials. 

Of course one could try to combine these two ideas. Thus one can compensate for a 
cutoff function R which depends on n by choosing the scale M to be /i. dependent too, 
chosen so that these spurious n dependent terms from M/ A modes in any physical result 
were removed. However it is extremely difficult to do this and M would have to be a 
complicated function of /x. In fact the only easy way to distinguish physical terms from 
spurious ones coming from M/ A modes is to compare against a calculation done with a 
regulating function R and renormalisation scale M, both of which chosen constant. In 
this case though it is clear that the one may as well just work with the fixed regulator 
prescription in the first place. 

Armed with these simple examples in terms of the cutoff A = M/A^/^, let me now 
return to the main calculations of this paper. The discussion in the previous section 4 
showed how to use the renormalisation scale M to relate all the different results but in 
several cases M had to depend on the physical parameter /i. This suggests that what 
is happening in the different calculations is that the UV modes around the scale M/r 
are being removed in different ways for different physical parameter values. Drawing on 
the comments in the appendix A. 2, it is clear that dimensional reduction is not the 
offender. Roughly speaking, dimensional reduction regulates by inserting a function 
Re — (/i;^/M^)~% which is completely independent of any physical parameters. 

However, ^-function regularisation regulates by altering the integrands. For instance 
for the L factorization, the regulator is, crudely speaking, J?^ = [{kf + (cu ± /i)^)/M^]~* 
in (3.10). This is clearly changing the way the UV modes are cutoff in a manner which 
depends on /j, and m. It is therefore no wonder that strange and m dependent terms, 
Y, appear even in the physical results. Equally, it should be no surprise that by choosing 
M to be a suitable, but complicated, function of /i and m, such terms can be removed. 

The lesson from the A cutoff examples is that one should not trust dependence on 
physical parameters in calculations using regulators which involve those same physical 
parameters. Thus in ^-function regularisation one should not trust the iioim dependence 
of the results for the free energy of a free Bose gas. Also, the only way to identify the true 
physics is to use a regulation scheme which is not sensitive to the physical parameters. 
Therefore one can not ascribe some fundamental meaning to one or other of the ^-function 
regularised forms results {pace [12, 13]). In the case of the free Bose gas one may as 
well just work with a simple constant UV cutoff or dimensional reduction as ^-function 
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regularisation merely adds unnecessary complications. 

6 Conclusions 

My first conclusion is to confirm the assertion of [9] that if and only if one includes 
Multiplicative Anomalies in ^-function regularisation calculations does one get mathe- 
matical consistency in the free Bose gas model. That is if one considers a well defined 
mathematical object, such as the (^-function regulated expression for F^^ (3.8), it does 
not matter then how you calculate it, either directly or via other quadratic expressions 
with the appropriate Multiplicative Anomalies, the same answer is always obtained. This 
is in contrast to suggestions elsewhere [13] that different calculational approaches to one 
well defined mathematical object in ^-function regularisation might give different answers 
because of some subtle physical effect. This point has been confirmed in several other 
models, see [2, 7, 14, 15] and appendix B for more details. 

The second conclusion of this work is that the ^-function regularisation calculations 
of the free Bose gas can not be trusted. Perhaps after a careful study of the regulation of 
the UV modes, one might be able to be sure that some (^-function regularisation results 
are free of spurious finite UV mode contributions. However it is far easier to use an 
physics independent UV regulation scheme, such as dimensional reduction, to find the 
correct physics in this case. Thus the thermodynamics of the free Bose gas is described 
by //3(m, T, //) of (3.7) alone, which is the standard result of the hterature e.g. [16]. The 
additional terms, such as Y{m, ji) of (3.5), found in some results, such as -F4^^, are genuine 
enough mathematically but they do not have a simple physical interpretation, they are 
merely regulation artifacts contrary to suggestions made in [9, 10]. 

This second conclusion means that one should not worry about the physics obtained 
when such artifacts are included in the free energy, as was done in [9] . All the different 
results, i^4i4'^, 6tc., all encode the physics of this model exactly but they do so in 
an extremely convoluted manner. Thus I also disagree with [12, 13] that only some of 
these expressions are physical, e.g. ^4/^:^, and the others, such as ^4^^, are unphysical. 
The dimensional regularisation calculations confirm that all the formal starting points are 
equally good. It is merely that some regularisation schemes are less convenient than oth- 
ers. Thus only by comparing ^-function regularised calculations against results obtained 
in good regularisation schemes, such as dimensional reduction, can one see that F^^kc, is 
the most convenient, but not more physical, starting point in ^-function regularisation. 

Thirdly, as was noted at the start, the extra terms in the free Bose gas results can 
sometimes be described in terms of Multiplicative Anomalies. Thus the results here throw 
some light on the importance of Multiplicative Anomalies in general. My results for this 
model confirm what has been noted elsewhere, namely that it is essential for mathe- 
matical consistency that Multiplicative Anomalies are included when using ^-function 
regularisation. However, the interpretation given here, and in the model used in [7], 
suggest that Multiphcative Anomahes have no novel physical content. They are merely 
contributions from high energy modes refiecting the way that ^-function regularisation 
cuts off these modes in a physical parameter dependent manner. 

Finally, these conclusions lead me to question the extraction of physical information 
from any ^-f unction regularisation calculation, even after Multiplicative Anomalies are 
accounted for. The analysis here was for one of the simplest examples, a free complex 
scalar field on a simple fiat Euchdean x space-time^^. Yet I have shown that it is 

^^The split of (3.2) and the role of UV divergences in this problem suggests that only the short distance 
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almost impossible to identify the correct physics from the (^-function regularisation results 
alone. It required a comparison with dimensional regularisation to do that. However, 
(■-function regularisation calculations are often used in much more complicated situations 
where surely it will be even harder to see the true physics from the UV regularisation 
remnants. I therefore believe that it will be almost impossible to identify the true physics 
rather than some regularisation remnants in many calculations using (-function regular- 
isation. 
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A Calculational methods. 

A.l Finite temperature calculations. 

First note that doing the Euclidean energy sum can be done using contour integration 
methods [18]. Suppose I have the following. 




(A.l) 



(A.2) 



The function g varies in this paper but one example is 




(A.3) 



no 



In this form, all the jj, dependence is in the propagators and hence in the gf's, and the 
energy sums are over pure Euclidean discrete values [19]. I use a subscript /? on integrals 
as a short hand for the Matsubara sum 

/ dh 9{h) = ^y2 9{kA = 27rn/3-^) (A.4) 

The contour C in (A. 2) is made of several pieces, each a small circle centred on the 
Matsubara frequencies z = 27rin/f3, running in the positive direction [18]. The one 
complication in the types of integral I am looking at is that there are cuts which may run 
along the real z (Minkowski energy) axis as the integrands are logarithms or polynomials 
in energy raised to some non-integer power. At the same time for the Bosonic fields being 
considered I have a Matsubara frequency at zero energy and therefore lying on the real 
axis. Thus I must distort the cuts to run either slightly above or slightly below the real z 
axis at z = 0. Normally this is trivial but the presence of a poles in the coth function at 
that point complicates matters there. Luckily, one can prove that bosonic functions, such 
as those relevant to the free energy expressions I am studying, have zero discontinuity at 
zero energy [24] so the cuts may be distorted to either side without changing the final 
answer. 

Now having dealt with that technicality, I distort the contour, expanding the contours 
round each of the poles on the imaginary energy axis until they merge and I am left with 
a contour running up the right hand side, and down the left hand side of the imaginary 
energy axis, i.e. from —too + e to from +ioo + e and then from +ioo — e to —ioo — e 
(0 < e <^ 1). Now I add on semicircular contours running at |2;| = cxd in each of the 
Re{z) > 0, Re{z) < half-planes. The integrand should be zero on these contours. The 
coth and tanh functions are constant there and usually the form I am given for f{z) dies 
away at large z. I am then sure that I have added nothing to the integral by adding these 
semicircular pieces to the contour. 

To make the result clearer I can rewrite the coth functions using 

1 1 ^Pz, 11 11 

-coth(— ) = +- + —5 = ^ , (A.5) 

2 ^ 2 ^ 2 e+^^ - 1 2 e-l^^ - 1 ^ ^ 

using the first (second) expression for the Re{z) > (<)0 part of the curve. This leaves 

/+fioo _ • 
—dz g{-iz) 
-ioo 27r 

F, = ^(-l)-^-(-°)Res[( ^^^J^^^_^ ,(-i.)| 



+ [ (-l)^Sn(.°)DisC 



1 



g/32*sgn(^o) _ \ 



(A.6) 



where {^o} are the set of poles of g. Res[. . .] and Disc[. . .] indicate that the residues at 
the poles and the discontinuities along the real z cuts should be calculated. The sgn 
function is given by sgn(2;) = -|-1 (—1) if Re{z) > 1 (< 1). 

The first term Fq is the Euclidean zero temperature energy integral. It has no explicit 
temperature dependence, though g contains chemical potential terms coming from the 
propagators. This shows that the UV finite T > contribution can be pulled off. All 
the discussion in this paper centres around terms which are independent of temperature. 
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though not of /i or m}^ Thus one can just focus on the UV infinite, T = terms 
Fq where one merely replaces the Matsubara energy sums in the original expression by 
a straight Euchdean integral.^° The chemical potential remains encoded here in the 
quadratic operators and propagators, present even in T = Fq expressions. 

One point of warning. A free Bose gas with non-zero charge density always has \fi\ = m 
at T = if Q/V 7^ since all particles must lie in the ground state in this case [16]. One 
must always remember that /j, is an implicit function of temperature since it is the charge 
density Q/V and not /i which is the physical observable. For a free Bose gas it may not 
always be possible to ensure < m for any given temperature and charge density. 

A. 2 Dimensional Regularisation 

In this appendix, I will summarize the most relevant aspects of the literature and fix 
my notation. References such as Collins [23] provide a much more detailed and careful 
derivation of the necessary mathematical theorems. However, 1 will also comment on how 
to apply this at finite temperature there is a crucial point which has direct relevance to 
the discussion of Multiplicative Anomalies. 

In dimensional regulation, one adds extra components to the loop momenta vectors in 
'directions' orthogonal to all other four-vectors. The number of these 'extra' dimensions 
is then taken to be — 2e, close to, but not exactly zero. More precisely the calculations 
are performed at a value of e where the integrals are well defined and then analytically 
continued to small finite e. Luckily there is no need to do these stages in detail as provided 
one follows some basic rules, dimensional regularisation can be implemented using some 
standard identities which guarantee mathematical consistency, at least for the simple field 
theories considered here. 

Suppose one has a single rf-dimensional Euclidean integration over a Euclidean loop 
momentum variable k with no external four vectors in the problem. If working at zero 
temperature and if the integrand depends only on K — k^, then dimensional regularisa- 
tion is implemented as follows 



where M is the renormalisation scale. Crudely, dimensional regulation controls the UV 
by inserting a regulating function R{k) of the form = (A;^/M^)^^, where e ^ Z, into 
the integrand. The UV divergences appear as 1/e terms. Note that IR divergences will 
appear in the same terms, a draw back for dimensional reduction in general but not an 
issue in the free Bose gas model. 

One useful identity in dimensional regularisation is 



^^I believe that one could easily make temperature dependent Multiplicative Anomalies given the 
explanation for Multiplicative Anomalies given elsewhere in this paper. All one needs to do is to make 
the UV regulating function temperature dependent. 

am assuming that one has kept the chemical potential in the propagators and not put it in the 
boundary conditions which is an alternative approach [18, 19]. 




(A.7) 





Demanding that this identity, and others hke it, are satisfied, fixes the forms of c{d) and 
g{e) [23]. For this to hold it is important that c{d) has the form given in (A. 7). 

For more comphcated integrals, e.g. with external loop momenta, consistency in di- 
mensional reduction requires that one regulates only using components of the loop mo- 
menta orthogonal to any external four- vectors in the problem [23]. When working in a 
heat bath though, there is the velocity of the heat bath with respect to the observer, 
u'^. I have chosen to work in the rest frame of the heat bath where u'^ = (1,0). Thus 

= k.u and, by the rules of dimensional regularisation, k^ can not be included in the 
regulating function. So I have in four-dimensions in thermal field theory for integrals 
with no external momenta and only one loop variable 

r d^k Dim.Reg. f d^'^'k M'^' _ f dkj f d^'H 

= c(3 - 2e)M2^i ^ j dK K^'^-\ (A.9) 
ki ■= neZ (A.IO) 

and c{d) is given in (A. 7). One can perform the energy sum in the usual manner assured 
that the result is finite even after the spatial integration provided one keeps e non-integer. 
One can then separate off the zero temperature and finite temperature parts as discussed 
in appendix A.l and shown in (A. 6). The former can then be returned to the usual zero- 
temperature (4 — 2e)-dimensional form through the use of (A. 8). The finite temperature 
part can have the regulator removed, e = 0, since it is UV finite. 

It is absolutely crucial that no regulation occurs in the energy variable k^ — k.u. So 
while the replacement 




{k^ = 2nn/P) (A.ll) 



may also regulate the divergent integrals, the resulting expressions will not in general be 
mathematically consistent. In other words, this type of regulation may have Multiplica- 
tive Anomalies! 

Finally, note that there is a close relationship between the family of different regulari- 
sation schemes and the families of renormalisation schemes. In dimensional reduction we 
can always change the definition of regularised integrals by multiplying by g{e), where 
g{z) is a function which is unit valued and analytic at 2; = 0. However, if one was to 
implement MS renormalisation but with a g dependent regularisation, i.e. just subtract 
the leading 1/e poles whatever g is, we would get extra finite terms proportional to g'{0). 
If is a simple algebraic function, and has no additional dependence on physical vari- 
ables, then these can be absorbed in constant redefinitions of the renormalisation scale 
M. For instance, in dimensional reduction it is common to use MS (Modified Minimal 
Subtraction scheme) which is equivalent to MS with a g{z) = 1 + 2;(ln(47r) — 7) + 0{z'^) 
and the relationship between these two schemes is simple. Alternatively, the effects of 
g{e) can be encoded through a conformal transformation on the dimensionless parameter 
e. 



A. 3 (^-function regularisation 



I implement ^-function regularisation scheme [4, 8, 25] as follows 

where a G Z+ is both the order and dimension of the operator A, i.e. I assume that 
A ~ I A; I" in the UV limit. The function g is set to one in standard calculations and 1 
will do that here, though see below for further comments about this. Roughly speaking, 
regulation is achieved by 



where = [A/Af J^^*'/*^. Thus both dimensional reduction and (^-function regularisation 
use non-integer powers, e or s, of k to achieve their regulation. The key difference for this 
work is the fact that ^-function regularisation's regulating function involves the operator 
A, and consequently varies with changes in any physical parameters in A. The regulat- 
ing function of dimensional reduction does not have this dependence. This ensures that 
dimensional regularisation satisfies certain basic identities [23], while (^-function regular- 
isation does not as appendix B notes. In fact this is the essence of the Multiplicative 
Anomalies in ^-f unction regularisation, while a great deal of work has been done to en- 
sure mathematical consistency of basic algebraic identities such as (A. 8) in dimensional 
regularisation [23]. A few further points are worth making. 

Firstly 1 have ensured that naively the expressions always have an overall factor of 
(M^)~* to mimic the (M^)~^ factor in the dimensional regularisation expressions. This 
aids comparison of the results obtained using the two regularisation schemes. More im- 
portantly though, it means that I have made a very specific choice for the form of the 
regulating power, which is not simply a small parameter s but a specific multiple of it. 
This is to ensure that the poles will cancel when calculating different terms in Multiplica- 
tive Anomalies as shown in (A. 15), i.e. I can use the same s parameter in all cases. In 
the same way, when comparing ^-function regularisation and dimensional regularisation 
expressions, a simple equality s — e will be sufficient to match such expressions. This 
is to be contrasted with the renomalisation scale M which will often be rescaled when 
comparing different expressions. 

Secondly, I choose to mimic dimensional regularisation in another way by ensuring 
that the UV divergence always appears as a s^^ pole with the physics contained in 
the 0(s°) term of a small s expansion. The generalised C-function C(-2|^) (for relevant 
operators A) is finite at ^ = but in general has poles at z = n + s,n = 1, 2, .... Thus 
creating pole is only an issue for the regularisation of logarithm of operators where 
I have had to add an overall factor of ^ in (A. 12). In fact this is very natural as the 
following example will show. 

Consider 

H = \n{AB) - \n{A) - \n{B) = (A.15) 

where A and B could be ordinary numbers or more complicated objects. These logarithms 
can be represented as the second term in the following series 

[{Ay/^] = 1 + ^ ln{A) + ^(ln(A))^ + 0{s') (A.16) 



The parameter a has no special meaning here but one can add it in anticipation of its 
role in the field theory case where A ~ \k"'\ as \k\ — >■ oo. Therefore I can construct a 
representation of H through the replacement 

ln(A) ^ -[{Ay/"] (A.17) 
s 

where s is small. Similarly for ln{B) I use s/b rather than s/a where b is arbitrary in this 
simple context. I must then use 

\n{AB) ^±^[(AS)*/('^+'')] (A.18) 
s 

to ensure that the l/s infinities cancel and I indeed find that 



2ab{a + b) 

which is zero as s — > 0. The first 0{l/s) terms in the expression cancel because of 
the judicious choice of regularisation parameters for (A.17) and (A.18). The next 0(s°) 
terms are just the original simple logarithms which, if we have constructed a sensible 
representation of the original expression (A. 15), must cancel. Thus the first non-zero 
terms are the 0{s^) terms but, provided all the terms are finite, this term will be zero in 
the limit s — where we expect to recover the expression H. However, if one is taking 
a trace over some variable upon which A and B depend, one can imagine that the trace 
may be divergent even if A and B are not. This may correspond to an additional l/s 
pole being generated which can then combine with these 0{s^) terms to give a non-zero 
contribution in the s — > limit, i.e. an anomalous H ^ result. Note the close similarity 
between the form of the expression for Hs and the Wodzkicki formula (A. 24) for the 
logarithmic Multiplicative Anomaly anomaly. 

One often hears that one of the advantages of (^-function regularisation is that it is 
finite. However all that is meant is that (^(2;|A) is finite near ^ = so one may define the 
Trlnyl = lndet74 = (^'{z = 0\A). In the context, of QFT this is not especially amazing. 
The physics is still in the second term of some expansion and one is merely giving a 
prescription for extracting this term. In my implementation of C-function regularisation, 
I have merely multiplied by a factor of l/s when regulating Trln expressions in order 
to keep other (^-function regularisation expressions and indeed dimensional reduction 
expressions on a similar footing. Thus the usual definition of ^-function regularisation, 
as applied to vacuum energy densities, corresponds to dropping what are the leading 
l/s poles in my expressions, i.e. traditional ^-function regularisation is also a 'minimal 
subtraction' renormalisation scheme. Ignoring the first term in the small s expansion 
is merely a particular renormalisation scheme, not just a simple regularisation^^. Put 
another way one can mimic the usual implementation of (^-function regularisation for 
vacuum energy densities in dimensional regularisation. To do this the usual dimensional 
regularisation expression is multiplied by e, which makes the expression finite in the same 
way that C(0|^) is. One would then define vacuum energies Trln74 = Indet A to be the e 
derivative of that combination, i.e. I take just the second term in a small epsilon expansion 

^"'^ Regularisation is process where infinities arc turned into large but finite terms. Renormalisation 
is a process which expresses bare parameters in terms of finite values which can be related to physical 
measurements. The latter must always be performed in QFT, even if the theory is finite, as quantum 
fluctuations are always present. It will though remove any infinities present in expressions for physical 
quantities. 
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+ 0{s') 



(A.19) 
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and drop the first term. This is of course just the Minimal Subtraction renormalisation 
scheme. 

Finally, one may always multiply the ^-function regularisation expressions, (A. 12) 
and (A. 13), by any function g{2s/a) provided this function g{z) is analytic and of value 
one at z = 0. Its introduction merely moves us through the family Schwinger proper- 
time regularisations [8], changing the finite terms by an amount proportional to g'{0). 
The g{0) = 1 condition can be imposed with out loss of generality as any other value 
is equivalent to a rescaling of the small regulating s parameter. This freedom is used 
to ensure that poles cancel in Multiplicative Anomalies, as discussed for equation (A. 17) 
and (A. 18). One example of interest is ^ = g{2s/a) — sT{s) for Trln problems so that 
from (A. 12) 

TrlnA'^-^:!^^-^c(^l^)r(.), (A.20) 

I will refer to this as SPTDR (Schwinger Proper Time Dimensional Regularisation). 
This member of the family of Schwinger proper time regularisation schemes is sometimes 
referred to as plain 'dimensional regularisation' in some of the literature on ^-function 
methods, as in [11]. However SPTDR is not the dimensional regularisation of the particle 
physics literature such as [17, 23] and as described in appendix A. 2. The review of Ball 
[8] makes this distinction clear and simple examples confirm this view. 

A.4 The Wodzicki residue. 

In the theory of elliptic pseudo-differential operators, there is a unique extension of the 
Diximier trace to elliptic pseudo-differential operators and Wodzicki gave the explicit 
form [5, 26]. This beautiful result of modern mathematics is central to several areas, 
such as Non- Commutative Geometry (e.g. see [6]). However, there is a simple formula 
for calculating it and a simple relation between the Multiplicative Anomaly of (1.1) and 
the Wodzicki residue. 

The simple recipe to find the Wodzicki residue^^ is as follows [5, 26, 1]. The "complete 
symbol" for our operators is A{x, k) :— e."^^^ Ae^^^ . Extract the UV behaviour as an 
asymptotic expansion 

— oo 

A{x, tk) = J2 t^M^^ k) (A.21) 

j=a 

where a is the order of the operator. The Wodzicki residue W of an operator A is then 

resw{A)= [ d^x [ 4^kM^S{k^ - M^)A^D{x,k) (A.22) 
Jm J 

Note that in many places (all the ones I found) the regularisation scale M is set to 
one, that is the calculations are done with all dimensionful parameters being measured 
in units of M. While perfectly acceptable mathematically, this does mean that the 
renomalisation scale is hidden even though it is extremely important when extracting 
real physical numbers from the mathematical results [17]. 

There is a simple link to the Cauchy residues of ^-functions through 

Res,=<,[C(s|A)] = a-heswiA-") (A.23) 

^^In [5] see sec. 7.13, pp.176. In [1] see section 6.5, pp. 225-226, especially his eqn.(9). Also note 
the comments in section 1.2, eqn.(4) which cites Wodzicki's 1984 thesis for a residue formula. See also 
Wodzicki [26]. 
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and finally to the Multiplicative Anomaly through 



a^{A, B) = {2ah{a + h))-h^w [{\n{A^ B-'')Y] (A.24) 

Note the strange expressions in (A.24) arc easily obtained in a simple analysis of loga- 
rithms of numbers which gives some insight into this expression, see (A. 19). 

A. 4.1 Conjecture for a generalised Multiplicative Anomaly formula 

The Wodzicki residue formula (A.24) is for Trln expressions (1.1), relevant to zero point 
energy calculations in general and to all the four-dimensional examples considered here. 
However simple traces over other operators can give rise to Multiplicative Anomalies, 
as the examples (B.2), (B.5) and (B.6) show. One can try to write these examples in 
terms of ii{z} = trln{y4}, with A — expj^;}, and then try to use the Wodzicki residue 
formula for the Multiplicative Anomaly. Unfortunately, the UV behaviour of the resulting 
operators A is then exponentially divergent or suppressed and the formula does not work. 
However, since direct calculations in (^-function regularisation in these cases does show 
the existence of an Multiplicative Anomaly, one wonders if there is not a formula similar 
to the one based on the Wodzicki residue. I have made the following conjecture. Suppose 
a Multiplicative Anomaly is given by an unregulated expression T 

•^ = E / Ai^. W, (A.25) 

where 

J2Fj{k) = V|A;| < oo (A.26) 
j 

The operators Fj in each term must have the same dimension which, without loss of 
generahty, I will take to be the same as the order a, that isF~|A;|"as/c— >oo. Let me 
now write down the ^-function regularisation form 

•^C = I] / ^''k [F,(A;)]^+^/"M-^ (A.27) 

where for dcfinitcness I choose to regulate the expression so that the UV behaviour is 
|^|o+s gj^^ ^YiQ renormalisation scale always appears as the simple M~* factor. Now 
suppose I make a double expansion of F^^^/"- in s and 1/k and define 



oo oo ^ 



m\n\ 

m=0 n=0 



The interest is in the UV divergences so it makes sense to separate off the integral over 
the I A; I through 

1 = / dt2t5{f -k^/W). (A.29) 
Jo 

By looking at the formula for the Tr In examples and trying a few examples, 1 have made 
the following conjecture 



oc 

j 



Y.'^jfk 5{l - P/M^)M».^^l-^/,,2,zp+„(^). (A.30) 



i.e. only the O(s^) and term of the fj expansion contributes. The reason for 

trying the n = D + a term comes from the original Multiplicative Anomaly formula using 
Wodzicki residue. There the term which appears is the one in the asymptotic expansion 
of the operator which is right on the boundary between divergence and convergence, i.e. 
the logarithmic divergent term which in this case is the n = D + a term. 

The reason for trying m = 2 term in the conjectured form (A. 30) is that it is quickly 
apparent that the m = and m = 1 terms are zero. That this must be so follows 
from the two constraints on Multiplicative Anomaly expressions. Firstly it should not 
itself be divergent so 1/s terms should cancel when doing the sum over j. In general 
one must choose the s parameter used to regulate each term in the j sum in just the 
right way to achieve this, and the definition used here ensures this, as noted in appendix 
A. 3. Secondly, the Multiplicative Anomaly expressions are naively zero in the first place 
(A. 26), and this is linked to the fact that the sum over j of the 0(s°) terms must be 
zero. Thus it can only be the third term in the s series of / that is relevant, presumably 
because in a proper calculation it is mixing with divergent terms. See (A. 15) for a simple 
example of this. 

B A collection of Multiplicative Anomalies 

In [7] it was suggested that Multiplicative Anomalies appear only if three criteria hold, 
namely 

1. the terms in an Multiplicative Anomaly must contain infinities, 

2. the regularisation scheme used must mix physical parameters (mass, chemical po- 
tential, etc.) in the artificial function used to cutoff UV momenta, and 

3. the different terms in the Multiplicative Anomaly must contain different combina- 
tions of these physical parameters. 

Thus regularisation schemes such as dimensional regularisation, physical cutoffs (e.g. 
lattice) and ones typically used in particle physics never have Multiplicative Anomalies, 
as the examples in [7] show. 

However, (^-function regularisation and all other members of the family of Schwinger 
proper time regularisations [8] are plagued by Multiplicative Anomalies. In general with 
these regularisation schemes regulated expressions will fail to obey many algebraic iden- 
tities naively satisfied by their unregulated counterparts. Is there any identity which 
C-function regularisation respects? 

The criteria above are merely necessary, not sufficient. For instance the simplest 
example of an Multiplicative Anomaly in QFT found in [2] is non-zero only in even 
dimensions from four upwards. However, numerous examples of non-zero Multiplicative 
Anomalies are now known. As some are relevant for the discussion of the free Bose gas, 
I will give a brief list. Here A, B are suitable operators, a is a c-number. 

The original example 

Tr\n{AB) ^ Trln(A) + Trln(5) (B.l) 

was given by Elizalde, Vanzo and Zerbini for several simple cases in field theory, including 
two real scalar fields of different masses in fiat Euclidean space-time [2] . Physically it is 
important when considering the vacuum energy of two free scalar fields as in cosmology 
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or in the Casimir effect. It also appears in the free Bose gas case, in particular the 
four-dimensional expressions (2.13), (2.19) and (2.22), which were first given in [9]. 

The problem of Multiplicative Anomalies appears with Green functions as well as 
with vacuum energy densities considered in [2]. Consider 

It is formally related to the first Multiplicative Anomaly, as defined in (1.1), by the re- 
placement in (1.1) A — > exp{l/y4} and B exp{l/{A+a)}. Examples of this failing with 
(■-function regularisation were shown in [7] for two real scalar fields of different masses^^ 
with A — A~^,a — 5m^. The Multiplicative Anomaly is then related to expectation 
values of the squares of the two real fields, {{(pf — (j)^) ^ — {(f^. 

Mass shifts are perhaps the most simple interaction possible and often appear in 
QFT. They are also afflicted by an Multiplicative Anomaly in (-function regularisation, 
for instance 

Tr{ln(A + a)] ^ Tr{ln(A)} - ^ iZ^Tr{A-"}. (B.3) 

n=l 

An example of this Multiphcative Anomaly was given in [7] for a single real scalar field, 
where one looks at a mass shift of the scalar field, A = A~^, a — 5w? 
A generalisation of the Multiplicative Anomaly in (B.l) is 

Tr{ln(n^n)} ^ J]Tr{ln(A„)}. (B.4) 

n n 

This version is relevant in the context of dimensional reduction [14, 15, 11]. 
Simple shifts by c-numbers lead to Multiphcative Anomahes 

TY{A + a} + Ty{A- a] ^Tx{2A}, (B.5) 

showing that the linear property of the trace is not preserved in (-function regularisation 
pace [20]. The first example of a Multiplicative Anomaly of this type was first noted by 
McKenzie-Smith and Toms [12, 13] in the context of the free Bosc gas using the three- 
dimensional expressions of (2.37) and (2.39). Again if in (1.1) I make the replacements 
A — > exp{— ^4 — a\ and B — > exp{— ^4 -|- a} this links this Multiplicative Anomaly to the 
original one. 

Finally, if one needed any further convincing about the existence of Multiplicative 
Anomalies, one need look no further than simple c-number multiplication of operators as 
in (-function regularisation 

Tr{a^} ^ aTx{A}. (B.6) 

It is needed as well as (B.5) when relating F3 and It is easy to see that if Tr{A} is 
regularised using C-function regularisation and becomes Tr{>l^+*} = a_is~^ + 00 + 0{s), 
then there is a Multiphcative Anomaly of 

arescaie(-4, o) := TrjaA} - aTr{A} = a\n{a)a-i. (B.7) 

This is non-zero whenever Tr{A} is infinite, just as required by the criteria in [7]. This 
is no more than a changing of the regularisation scale which always appears in QFT e.g. 
let A = A the propagator for a real scalar field, one will have to remember that I need 
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The analysis in [9] confirms that of [7] contrary to what is stated in [20] 



to introduce a new scale, say M, to keep the dimensions of the regularised expression 
correct, e.g. as in (3.23) and (3.9). Note that this tells us that 'trivial' operations such as 
pulling out the factor of two on the right hand side of (B.5) or pulling a minus sign outside 
the trace are in fact highly non trivial if Multiplicative Anomalies are non-zero. Yet such 
operations are performed frequently in formal derivations, and in particular they abound 
in our canonical derivations of the free energy of the free Bose gas in the last section. 
Thus the Multiplicative Anomaly, agMft of (2.40), appearing used in the three-dimensional 
forms for the free energy encodes a failure in a combination of (B.5) and in (B.6). 

The point of giving this summary has been to show how many different algebraic op- 
erations fail when using ^-function regularisation, even in the simplest of QFT problems 
involving free fields in fiat Euclidean space times. One can make these algebraic opera- 
tions but only if one remembers to include the Multiplicative Anomaly terms as well. It 
is extremely easy to neglect these Multiplicative Anomahes but this will lead to mathe- 
matical inconsistencies. In this sense Multiplicative Anomalies are of vital importance to 
problems using ^-function regularisation. 

Several of these Multiplicative Anomalies are relevant to the analysis of the free energy 
of the free Bose gas. In particular, the regularisation scale, here called M, is often shifted 
in the analysis above. The simple example (B.7) can be interpreted in this way. This 
all gives further weight to the idea in [7] that one can take account of all Multiplicative 
Anomalies by regarding them as shifts in renormalisation scales. 

C Useful formulae 



n=0 
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Jo 



oo x^~^ 



dx- — — = , (C.2) 

{i + i3xpy p r(i/) ^ ^ 



r(n + e) = r(n) (1 + eV'(n)) , n e Z+ (C.3) 
1 (—1)"- 

r{-n + e) = -J-J— (l + ei;{l + n)), n G {0, Z+} (C.4) 



r(l/2) = ^Ar (C.5) 
2 4' 



r(3/2) = Wn)-^ (C-6) 



r(-l/2) = -2T{1) = -2xAF (C.7) 



tlj{l + x) = tlj{x) + - (C.9) 

X 

^(1) = -7 « -0.57721566490, 7 is Euler's constant, (C.IO) 
^(1/2) = ^(l)-21n(2) (C.ll) 



00 



j^^K = ''^^27r dKK''-' £ de{sin{e)f-^ (C.12) 

POO 

= c{D) / dKK^-\ (C.13) 

= ft kd = Kcos{9), \k\ = Xsin(^) (C.14) 

c(d) = ((47r)'^/'r(d/2))"' (C.15) 
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